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INTRODUCTION  AMP  BACKGROUND 

Optimal  periodic  control  has  become  an  Increasingly  Important  area  of 
research  as  evidenced  by  the  greater  numbers  of  published  papers ,  theses, 
and  dissertations  resulting  from  government  grants,  academic  Interest  and 
Industrial  sponsorship.  Examples  of  periodic  or  cyclic  processes  abound 
In  nature,  ranging  from  the  rhythmic  pulse  o*  living  creatures  to  the 
perpetual  orbits  of  celestial  bodies.  Many  are  examples  of  an 
optimization  process  of  nature.  In  contrast,  there  are  numerous 
engineering  systems  designed  to  operate  in  a  steady  Btate  conditions.  In 
many  of  these  systems  performance  could  be  aignif icantly  improved  by  some 
form  of  cyclic  operation.  Chemical  plant  process  control  was  one  of  the 

first  to  be  investigated  for  improvement  compared  to  steady  state 

12  3 

performance  ’  *  .  This  motivated  the  first  paper  on  optimal  periodic 
4 

control  and  led  to1  its  subsequently  rapid  theoretical 
5  6  7  8 

development  '  '  ’  .  Sur-ey  papers  V,)0,  and  11  sumnurize  major  results 
through  1975. 

The  research  encompassed  by  thi9  report  was  stimulated  by  a 
controversy^'^’^  over  an  aerospace  problem.  The  nonoptimality  of 
sceady  state  cruise  for  an  aircraft  with  respect  to  fuel  efficiency  was 
shown  in  1976  by  Speyer^.  However,  subsequent  effort*^  to  find  a 
locally  optimising  solution  to  this  problem  failed  using  standard 
optimization  computational  techniques  such  as  steepest  ascent  and 
conjugate  gradient  methods.  This  led  Speyer  and  Evans  to  the  formulation 
of  a  minimum  state  optimal  periodic  control  problem  that  would  generate 
periodic  solutions.  A  locally  optimizing  solution  to  this  new  problem  was 


obtained  in  the  form  of  an  asymptotic  expansion  about  a  small  parameter. 


This  analytical  result  provided  excellent  agreement  with  numerical  results 

obtained  In  a  parallel  effort17.  Further  numerical  study  of  this 

problem  and  the  development  of  corrections  and  additions  to  the  general 

1 8 

theory  of  optimal  periodic  control  followed  .  The  results  of  this  work 
provided  a  better  understanding  of  the  source  of  failures  previously 
encountered  In  attempts  to  solve  the  aircraft  cruise  prohlem. 

This  report  describes  a  continuation  of  this  research  conducted  at  the 
Frank  J.  Seiler  Research  Laboratory  under  Work  Unit  2305-F2-67.  The 
principal  results  of  this  effort  have  been  presented  at  several 
conferences  and  have  bean  published  In  thq  open  literature.  The  major 
accomplishments  obtained  are  summarized  in  the  remainder  of  this  report 
under  three  headings:  (1)  A  Second  Variation  Condition;  (2) 

Computational  Techniques;  and  (3)  Optimal  Aircraft  Cruise. 

This  research  has  been  a  collaboration  of  work  by  the  orlnclpal 


Evans ,  Lt  Col,  USAF;  Professor  Jason  L,  Speyer. 


The  University  of  Texas  at  Austin;  and  two  of  his  graduate  students,  David 


E.  Walker  and  David  P.  Dannemlller.  The  work  by  Professor  Speyer  and  his 


students  was  partially  sponsored  by  the  National  Science  Foundation  Grant 


ECS7918246.  The  activities  of  the  students  also  contributed  significantly 


to  their  Master's  degrees. 

A  complementary  research  effort  was  sponsored  by  AFOSR  Grant  Number 

77-3158  during  the  period  1  October  lq76  to  31  January  1982.  Personnel 

associated  with  this  work  included:  Principal  Investigator,  Professor 

Elmer  G.  Gilbert,  The  University  of  Michigan;  Arthur  E.  Frazho, 

post-doctoral  researcher;  and  PhD  students  Daniel  T.  Lyons  and  Dennis  S. 

19 

Bernstein.  See  Professor  Gilbert's  Final  Report  for  a  summary  and  a 


complete  bibliography  of  their  work  In  optimal  periodic  control. 


OBJECTIVES  OF  RESEARCH 


"O)  The  slated  objective  of  this  research  was  to  develop  the  theory  and 
computational  technique  for  optimizing  the  flight  path  of  an  aircraft  with 
respect  to  fuel  consumption  ^maximize  range  for  a  given  amount  of  fuel) 
during  the  cruise  segment  of  flight.  It  then  was  Intended  to  apply  these 
tools  to  a  point  mass  model  of  an  aircraft  and  determine  a  locally 
optimizing  cyclic  cruise  flight  path  In  a  proof  of  principle  demonstration. 

Their  are  numerous  potential  applications  for  this  research.  The  more 
obvious  Air  Force  benefits  include:  extend  the  range  of  an  air  vehicle 
with  a  fixed  amount  of  fuel;  reduce  Its  fuel  requirements  for  a  given 
range  thereby  Increasing  Its  load  capability;  and  Increase  Its  endurance 
allowing  It  to  remain  aloft  longer.  The  Improvements  of  periodic  cruise 
flight  paths  appear  to  be  most  suited  for  remotely  piloted  aircraft  or 
cruise  missile  type  applications.  However,  in  many  emergency  or  back-up 
operations,  it  also  would  be  quite  feasible  for  manned  systems. 


PROBLEM  FORMULATION  AND  FIRST  ORDER  NECFSSARY  CONDITIONS 

The  optimal  periodic  control  problem  consists  of  minimizing  the 
performance  criterion 


J(u(»).  *0>  T  ) 


if' 

i  Jo 


L(x(t),u(t))dt, 


(1) 


with  respect  to  the  perlo •'  T£  T  =  (O,00),  the  p-vector  control  functions 
u(  •  )  ell,  where  U  Is  defined  in  Assumption  2  below,  and  the  initial  states 
x(0)  ■  Xq  e  Rn,  subject  to  the  time-invariant  dynamical  equations 


x(t)  ■  f (x(t) ,u(t)) , 


(2) 
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villi  the  periodic  boundary  conditions 


x(0)  -  x(T).  (3) 

Note  that  both  the  integrand  of  the  performance  index  and  the  dynamical 
equations  are  tirae~invariant . 

Assumption  1:  f(*,*)  and  !.(*,•)  and  their  derivatives  up  to  second 

order  are  assumed  to  be  continuous  with  respect  to  both  arguments. 

Assumption  2:  U  =  ( u (  ^  ):  u(  •  )  Is  piecewise  continuous  in  the 

interval  (0,°°)  and  ||u(»)||<*1  =  sup  |u(t)|  <®  where 

|u(t)|  =  (t  ,.^(t))1/2  ,  u(t)  eRpJ 

i“l 

Definition :  A  piecewise  continuous  function  f(*  )  has  a  period  if 

there  exists  a  minimum  te  T  auch  that 

f(T)  -  f(0)  (4) 


This  minimal  T  Is  caLled  the  period  of  f(»). 

Remark:  This  definition  excludes  constant  f(*). 

The  first  order  necessary  conditions  for  optimality  derived  from  the 
calculus  of  variations  are: 


f (x,u) , 


(5) 


(6) 
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(7) 


0  -  H 


x(0)  -  x(t). 


X(0)  -  A(x), 


H(t)  -  J(u(.).x0iT)  -  0, 


where  X(t)  t  Rn  is  the  Lagrange  multiplier  that  adjoins  the  system 
constraints  (2)  to  the  performance  index,  L(x,u),  forming  the  variational 
Hamiltonian  defined  as 


H(x , u, X)  ■  L(x,u)  +  f(x,u). 


It  Is  also  assumed  that  the  Legendre-Clebsch  condition  is  met  in  strong 
form  along  the  extermal  path,  i.e., 


H  >  0. 
uu 


Any  periodic  solution  to  the  two  point  boundary  value  problem, 
equations  (5)  through  (9)  is  an  extremum  of  the  problem.  The  condition 
(10)  relating  the  Hamiltonian  and  the  performance  Index,  evaluated  along 
the  optimal  path  Is  the  special  condition  for  testing  the  optimal  period, 

4 

firet  derived  by  Horn  and  Lin  . 
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A  SECOND  VARIATION  CONDITION 


Solutions  that  satisfy  the  first  order  necessary  conditions,  equations 
(5)  through  (10)  end  (12)  are  examined  for  local  optimality  by  second 
variation  tests,  ouch  as  the  Jacobi  necessary'  condition.  A  very  useful 

form  of  the  Jacobi  test,  developed  by  Bittanti  et.  al.^  and  extended  by 

20 

Gilbert  and  Bernstein  ,  shows  whether  or  not  a  static  solution  is 

locally  optimal.  This  test  was  used  by  Speyer^"*  and  by  Brcakwell  and 
21 

Shoaee  to  show  that  static  cruise  for  many  aircraft  models  is  not  fuel 

minimizing.  However,  this  test  provides  insufficient  information  to 

determine  the  optimality  of  cyclic  or  periodic  solutions. 

An  important  result  of  the  research  covered  hy  this  report  is  the 

development  of  a  variational  theory  for  testing  periodic  solutions.  This 

22 

work  was  ’--esented  by  Speyer  at  the  1981  Joint  Automatic  Control 

Conference  in  Charlottesville ,  VA  and  will  be  published  this  fall  in  the 

23 

IEEE  Transactions  on  Automal  lc  Control  .  Two  earlier  papers,  clarified 
and  extended  by  this  effort,  are  Brlltanti,  et.  al.^-  who  considered  the 
problem  for  fixed  period,  and  Chaivg^,  who  extended  the  work  to  free 
period.  The  new  results  are  summarized  in  the  remaining  paragraphs  of 
this  section. 

Properties  of  autonomous  Hamiltonian  systems  and  their  related 
monodromy  matrix  are  used  to  establish  relationships  essential  to 
developing  the  second  variational  conditions  for  optimality.  The 
monodromy  matrix  is  the  transition  matrix  for  the  state  equations  (5) 
through  (9)  evaluated  over  one  period.  It  is  shown  that  the  monodromy 
matrix  has  at  least  two  unity  eigenvalues  and  that  two  cf  them  are  coupled 
(in  the  same  Jordan  block).  The  eigenvector  associated  with  one  of  the 


unity  eigenvalue*  la  tangent  to  the  state  space  orbit  described  by  the 
related  ext renal  solution.  The  generalised  eigenvector  associated  with 
the  other  unity  eigenvalue  defines  the  direction  of  a  one~d laenslonal 
family  of  orbits  which  varies  continuously  with  the  Hamiltonian. 
Determining  this  direction  was  a  key  factor  in  developing  the  algorithm 
used  for  computing  eztremal  aolutlona  to  the  optical  periodic  control 
problem  described  in  the  next  section. 

Another  important  relationship  derived  from  this  effort  Involves 
classifying  the  eigenvalues  that  can  result  for  real  values  of  the  Riccatl 
variable.  The  existence  of  a  rsal-vaiued  solution  to  the  Riccatl 
differential  equation  Is  a  veil  known  second  variation  condition.  When 
the  solution  to  the  Riccatl  equation  is  periodic  it  can  be  expressed  In 
the  form 


P*l2P  +P*  11  -  *22P 


21 


-  0, 


(13) 


where  P  la  a  vector  of  Initial  conditions  defining  the  periodic  solution 
and  the  t'e  are  square  partitions  of  the  monodroay  matrix.  The  canonical 
similarity  transformation  of  the  monodromy  matrix  gives 


L  *L 


-1 


*11  +  *12P 
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*22  "  ?*12  J 


(14) 


where  the  Identity  (13)  is  used  to  obtain  the  rero  element  and  the 
similarity  transform  matrix  is  specified  as. 


7 


The  eigenvalues  of  this  transformed  matrix  (14)  are  those  of  the  monodromy 
matrix  due  to  their  Invariance  through  a  similarity  transformation. 
Because  of  the  sero  element  on  the  off-dlagonal  of  the  transformed  natrlx, 
the  eigenvalues  of  the  submatrices  on  the  diagonal  are  the  same  as  for  the 
entire  matrix. 

Osing  the  symplectlc  property  of  the  monodromy  matrix,  which  is  also 
preserved  through  the  similarity  transformation,  the  following  important 
relationship  is  obtained, 


The  significance  of  this  equation  is  that  it  strongly  restricts  the 
eigenvalues  of  the  monodromy  matrix  that  correspond  to  real-valued  Riccatl 
variable  elements  since  the  elements  of  the  monodromy  matrix  must  alec  be 
real-valued  for  a  physically  realisable  system.  Recall  also  the  matrix 
property  that  the  determinant  of  a  matrix  la  equal  to  the  product  of  its 
eigenvalues.  Considering  eigenvalues  of  magnitude  one,  the  following 
result  can  be  stated: 

A  necessary  condition  that  the  Riccatl  variable  matrix,  P,  be  real¬ 
valued  Is  that  there  be  no  distinct  eigenvalues  of  the  monodromy 
matrix  on  the  unit  circle. 
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Remark:  The  satisfaction  of  this  condition  does  not  guarantee  that  P 
exists  for  all  starting  times  over  the  Interval  of  a  period.  The  solution 
of  the  Rlccatl  differential  equation  is  still  required  to  ensure  that 
there  are  nc  finite  escape  Claes. 

A  second  variational  sufficiency  condition  for  weak  local  optinallty 
of  cyclic  processes  can  also  be  stated.  The  following  condition  extends 
and  clarifies  previous  statements^' ^  of  the  condition: 


For  the  periodic  control  problem  described  by  (1)  through  (3)  and 


assumptions  1  and  2,  (u°(  •  ),  x°,  f°)  £  U  x 

0 

local  ainlaum  if; 


x  T  forms  a  weak 


(i)  the  first  order  necessary  conditions  (5)  through  (10)  are 
satisfied, 

(ii)  the  strong  form  of  the  Legondre-Ciebsch  condition  (12)  is 
satisfied, 

(111)  there  exists  a  real  valued  bounded  symmetric  matrix 
solution  to  the  Riccati  differential  equation  on 
0  <  t  <  T°  satisfying  the  periodic  condition 
P(0)  -  P(ip), 

(iv)  there  are  no  eigenvalues  of  the  monodromy  matrix  on  the 
unit  circle  except  for  the  two  coupled  unity  (+1) 
eigenvalues  associated  with  the  velocity  vector  y(0) 
where  it  (H)/ffi  t  0  ensures  this  coupling,  and 

(v)  tne  eigenvalues  of  Che  mouodromy  matrix  off  the  unit 
circle  are  distinct. 
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Remarks : 


1.  The  requirement  for  earlier  statements  of  this  sufficiency 

condition’’’^  chat  the  matrix  ♦  (  T  ,0)=  ♦  have  no  unity 

eigenvalue  la  never  satisfied.  Furthermore,  it  la  required  in  condition 
(lv)  that  the  remaining  eigenvalues  of  *^(1,0)  not  be  on  the  unit  circle. 

2.  Condition  (v)  la  a  fora  of  the  strongly  positive  condition  which 

is  totally  lacking  in  the  previous  stateaents  of  this  sufficiency 

condition. 

COMPUTATIONAL  TECHNIQUES 

Although  periodic  optimal  control  problems  had  formed  an  Important 
class  of  practical  problems,  few  numerical  investigations  had  been 
reported  through  1980.  Initial  experimentation  indicated  poor  convergence 
behavior  for  first  order  optimisation  schemes  relative  to  a  class  of 

1  fi 

aircraft  cruise  problems.  It  has  been  recognized  that  this  was  due  in 
part  to  the  shallow  curvature  of  the  cost  criterion,  partly  due  to  the 
lack  of  eenaitivity  of  first  order  methods,  and  finally  due  to  the  great 
difficulty  in  closing  the  solutions  (satisfying  the  periodicity 
requirements) . 

In  order  to  develop  a  better  understanding  of  the  difficulties 
encountered  In  earlier  numerical  Investigations  of  the  optimal  cyclic 
aircraft  cruise  problem,  an  Illustrative,  alulaua  state,  optimal  periodic 
control  problem  was  formulated.  An  analytical  solution  to  this  problem 
was  first  obtained^  using  a  perturbation  mathod  most  frequently 
credited  to  Llndstedt  and  Poincare.  The  solution  can  be  expressed  in  the 
form  of  an  asymptotic  series  expansion  about  a  small  parameter  which  also 


10 


can  be  written  In  the  fora  of  a  Courier  aeries  expansion.  This  result 

captures  an  interesting  characteristic  of  the  solutions  that  satisfy  the 

first  order  necessary  conditions,  equations  (5)  through  (9).  That  is, 

they  fora  a  set  of  solutions  varying  continuously  in  amplitude  and 

period.  The  results  of  this  study  were  presented  by  Evans^  at  the  1979 

Joint  Automatic  Control  Conference  in  Denver,  CO.. 

18 

The  initial  numerical  study  of  the  illustrative  problem  showed 

that  the  approximate  analytical  solution  was  quite  good.  It  verified  the 

infinity  of  solutions  that  satisfy  equations  (5)  through  (9)  and  form  a 

continuous  set  or  family.  However,  it  was  also  discovered  that  an 

infinity  of  families  of  solutions  were  found  to  exist.  The  families 

Intersect  at  common  solutions  called  bifurcation  points.  An  illustration 

of  these  results  Is  given  later  in  the  section.  It  should  be  noted  that 

much  work  has  been  accomplished  by  dynamlclsts  determining  periodic 

solutions  to  a  set  of  first  order  differential  equations.  The  work  by 
24  25 

Henon  and  by  Contopoulos  has  been  invaluable  in  this  research  and 
it  provides  a  detailed  characterization  of  solution  families  and 
bifurcation  points. 

The  results  from  the  numerical  investigation  of  the  illustrative 

26 

problem  were  presented  by  Evans  at  Che  2nd  International  Federation  of 
Automatic  Control  Applications  of  Nonlinear  Programming  and  Optimization 
at  Oberpfaffenhofen,  West  Germany  in  September  1980.  The  emphasis  here 
was  placed  on  identifying  the  richness  and  complexity  of  the  solutions  to 
this  type  of  control  problem.  Characteristics  of  solutions  and  of 
families  of  solutions  were  examined  in  some  detail.  The  associated  paper 
was  published  in  the  Conference  Proceedings. 
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Ttys  computer  program  chat  was  used  for  Che  initial  study  was  an 

27 

adaption  of  on«  developed  by  Broucke  to  find  periodic  solutions  to  4th 

order  dynamic  systems.  Search  methods  about  a  known  solution  were  used  to 

find  new  solutions  of  the  family.  As  indicated  in  an  earlier  section  the 

direction  of  the  family  can  be  predicted  from  the  generalised  eigenvector 

associated  with  the  second  unity  eigenvalue  of  the  monodromy  matrix.  A 

shooting  method,  using  these  predicted  starting  values,  Integrates  the 

Euler-Lagrange  equations  (5)  through  (9)  to  obtain  periodic  solutions  in 

convergent  iterations.  The  improved  shooting  method  and  additional 

results  from  the  numerical  investigation  of  the  problem  was  presented  by 

2s 

Evans  at  the  20th  IEEE  Conference  on  Decision  and  Control  at  San  Diego 
in  December  1981.  The  associated  paper  was  published  in  the  proceedings 
to  the  conference. 

The  remainder  of  this  section  expresses  important  concepts  and  key 
relationships  associated  with  the  development  of  the  computational 
technique  employed  during  this  research  effort.  First,  identification  of 
solutions  is  most  easily  accomplished  by  association  with  their  initial 
conditions  since  a  set  of  initial  conditions  Identifies  a  unique 
solution.  An  lmportent  characteristic  of  a  periodic  solution  is  the 
number  of  axis  crossings  in  the  same  direction  that  occurs  during  one 
period  of  the  solution  for  a  particular  variable  of  the  problem.  This  is 
a  distinguishing  characteristic  of  a  family  of  solutions.  In  most  cases, 
the  number  of  axis  crossings  is  the  same  for  all  solutions  of  the  family. 
This  can  be  used  as  a  program  check  to  verify  that  a  new  solution  belongs 
to  the  family  that  was  Intended  to  be  followed. 
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After  one  solution  has  been  obtained,  the  initial  conditions  for  the 
next  solution  on  the  family  can  be  projoced.  With  the  definition, 

¥  =  y(0)  -  y(t )  (16) 

a  small  change  In  ¥  due  to  a  variation  in  the  Initial  conditions  y(0)  and 
the  period  (T)  gives  the  following  results: 

d  V  ■  Sy(0)  -  6  y(t)  ■*  y(x)  dy 

■"  “  ,f(T>  *'1'  [  «y<0>] 

Two  elements  must  be  fixed  to  use  to  predict  new  guesses  since  (1) 

y(T)  le  proportional  to  a  column  in  ^-I,  and  (2)  $-1  becomes  singular 

as  a  solution  is  approached.  Removing  corresponding  columne  of 
[y(t)  ^-1]  eliminates  lndetermlnancy  resulting  in 

-d  Hf  -  Sl6z  (1.9) 

Taking  the  pseudo  inverse  of  0  allows  computation  of  new  smarting 
conditions  from  {z  for  determining  a  new  solution  of  the  family.  A  less 
cumbersome  predictor,  such  as  a  curve  fitting  interpolator,  la  suggested 
after  several  solutions  have  been  obtained. 

To  Illustrate  the  computational  technique  some  results  of  the  minimum 
state,  optimal  periodic  control  problems  will  be  uaad.  The  Euler  Lagrange 
equations  derived  from  the  first  order  conditions  for  this  problem  are 


(17) 

(18) 
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(20) 


*1  *  *2  ' 

«  *2 
*2  "  '  T 

\  "  -*1  * 

1  3  X 

A2  -  *2  -  *2  l  ♦ 

and  x(0)  -  x(T)  ,  X(0)  -  X(  T) .  (24) 

For  this  problem  a  static  solution  exists;  l.e.,  x^  -  x2  “  X^  “  X  2  “  ® 
Starting  from  this  solution,  a  family  of  solutions  can  be  determined  as  In 
Figure  1.  The  solutions  are  represented  by  initial  conditions.  For  this 
illustration  x^O)  «  X^(0)  •  0  for  all  points  graphed.  The  Initial 

value  of  is  represented  on  the  graph.  The  last  CGuditior.  is 

determined  by  the  relationship  of  the  other  states  and  the  Hamiltonian 
evaluated  at  the  initial  time.  The  state  relationships  (x^  vs  x2)  for 
several  solutions  are  superimposed  on  the  graph  of  the  family  and  centered 
at  points  corresponding  to  their  respective  initial  conditions.  The  scale 
of  the  x^  vs.  x2  plots  are  all  the  same. 

As  indicated  before,  an  Infinity  of  families  of  solutions  exist. 

Plots  similar  to  the  previous  one  are  depicted  lc  Figures  2,3,  and  4  for 
three  additional  families.  Each  plot  is  to  the  same  scale  as  In  Figure 
1.  Note  the  number  of  axis  crossings  in  the  various  examples.  For  a 
single  family  the  number  cross Inga  are  generally  the  same.  The  family 
that  emanates  from  the  static  solution  is  called  the  principal  family;  all 


(21) 

(22) 

(23) 
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other*  ere  branch  families.  The  solution  identified  by  c  In  Figures  1 
through  4  represents  the  solution  of  the  respective  family  that  also 
satisfies  the  optimal  period  conditions,  equation  (10). 

Several  additional  levels  of  brai  'hlng  are  shown  In  the  detail  of 
Figure  5.  Note  that  families  branch  only  In  "stable’'  areas  of  the 
family.  Here  stability  refers  to  no  eigenvalues  of  the  monodromy  matrix 
existing  outside  the  unit  circle.  Bifurcation  points  (branch  points  of 
the  families)  are  dense  In  the  stable  regions  of  the  family. 


•  o  01  9.}  0.1  0, 


Figure  1.  Variation  of  Periodic  Solution  Along  Principal  Family 
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Figure  4.  Variation  of  Periodic  Solution  Along  Branch  Fa«lly  F5/2A 
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bifurcation  points. 

a.  2(S/ZA) 


B.  lO/ZA)  C.  2  (1  \5/ZA) ) 


Figure  5.  Branch  Families  (Detail) 


The  purpose  of  this  detsll  is  Intended  to  Identify  some 
charscterlstics  of  the  solutions  snd  to  emphasise  the  necessity  of 
developing  s  systematic  approach  for  Investigating  optimal  periodic 
control  problems.  The  amplitude  and  period  of  solutions  along  a  family 
vary  continuously.  However,  sharp  dlfferencea  in  the  amplitude  and  period 
of  solutions  of  different  families  generally  exist  even  though  Initial 
conditions  for  both  may  vary  only  slightly.  See  reference  18  for 
Identifying  bifurcation  points  and  branch  families. 

Suosaarlzlng  the  computational  technique;  first,  a  shooting  method  Is 
used  to  find  a  closed  periodic,  path  which  satisfies  the  first  order 
necessary  conditions  (the  Euler  Lagrange  equation  (5)  through  (9),  except 
for  the  tranaversallty  condition  (10)  associated  with  free  period).  Then 
a  one-dimenslonal  family  of  periodic  solutions  Is  constructed  using  the 
generalized  eigenvector  or  a  curve  fitting  Interpolator  to  predict  initial 
conditions  of  additional  solutions.  Finally,  the  family  is  traversed  In 
the  direction  of  decreasing  cost  criterion  until  the  optimal  period 
condition  (10)  Is  satisfied. 

OPTIMAL  AIRCRAFT  CRUISE 

Fuel  efficient  cruise  trajectories  for  aircraft  have  been  a  subject  of 
continuous  theoretical  Interest  and  are  becoming  one  of  practical  Interest 
as  well.  Since  the  steady  state  cruise  path  la  not  mtnlmitlng15  for 
most  point  mass  aircraft  models,  the  objective  Is  to  obtain  the  periodic 
paths  that  are  minimising.  There  appear  to  be  two  underlying  mechanisms 
for  producing  periodic  paths.  The  first  mechanism  is  the  mismatch  in  the 
regions  of  velocity  and  altitude  where  the  aircraft  is  aerodynamic  and 
propulsion  efficient.  This  Is  the  mechanism  behind  chattering  (or  relaxed 


steady  state)  cruise.  There  also  is  a  potential  and  kinetic  energy 
interchange  which  is  optimal  for  fuel  interchange  which  is  optimal  for 
fuel  performance.  The  need  for  substantial  kinetic  energy  seems  to  be  the 
reason  for  the  velocity  threshold  found  earlier15. 

29 

Recent  work  in  thic  area  includes  the  Master' o  Thesis  of  Walker 
and  Dannemiller30.  Both  applied  the  techniques  summarised  in  this 
report  to  investigate  the  optimal  periodic  cruise  of  a  hypersonic 

cruiser.  The  results  of  this  work  were  also  presented  by  Speyer31  at 
the  AIAA  Guidance  and  Control  Conference  at  Danvers,  MA,  in  August  1980. 

A  point  mass  model  of  an  atmospheric  vehicle  operating  in  the 

hypersonic  region  was  used  to  investigate  the  fuel  improvement  from  the 
steady  state  cruise  path  obtained  by  modulating  the  flight  path.  The  fuel 
improvement  obtained  was  due  solely  to  a  potential-kinetic  energy 
interchange  which  was  indicated  by  a  frequency  type  second  variational 

analysis  of  the  steady  state  cruise  for  the  flat  earth  model.  A  family  of 
solutions  was  generated  for  both  the  flat  and  spherical  earth  models.  By 
applying  the  second  variational  sufficiency  conditions  for  periodic 
processes,  only  one  flight  path  which  involves  the  flat  earth  model  was 
found  to  be  locally  minimizing.  The  improvement  of  the  periodic  cruise 
over  the  steady  state  cruise  for  this  example  is  4.2Z.  Ho  locally 

minimizing  path  for  the  spherical  earth  model  was  found.  Nevertheless, 
the  periodic  extremal  cruise  paths  found  did  improve  fuel  performance  over 
their  respective  steady  state  cruise  paths  by  as  much  as  4.5Z. 


RECOMMENDATIONS 


The  three  principal  objectives  of  this  research  have  been 
satisfactorily  achieved.  Even  though  the  minimum  state  Illustrative 
problem  has  been  exhaustively  studied,  further  Investigation  has  merit. 
Verification  of  new  theory  or  new  computational  techniques  are  more  eaoily 
accomplished  with  the  reduced  state  problem.  Additional  relationships  may 
be  exhibited  by  further  investigation  of  out-of-plane  solutions, 
eigenvector  directions  at  bifurcation  points,  and  other  solutions  that 
also  might  satisfy  all  first  and  second  order  conditions. 

Now  that  a  locally  optimising  periodic  flight  path  has  been  found  for 
one  model,  potential  applications  should  be  examinee.  More  realistic 
aircraft  models  should  be  developed  and  st  .idled.  Certainly,  the 
feasibility  of  subjecting  the  aircraft  to  cyclic  control  oust  be 
considered.  In  particular  the  cycling  on  and  off  of  its  engines. 

A  related  area  of  research  Is  associated  with  quasl-perlodlc  solutions 
which  may  provide  better  performance  In  some  Instances  than  t  e  periodic 
solutions.  A  considerable  amount  of  research  effort  In  this  area  has  been 
expended  by  statistical  dynaolclsta. 
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